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In this report, we discuss a candidate mechanism through which one might address the various 
cosmological constant problems. We first observe that the renormalization of gravitational couplings 
(induced by integrating out various matter fields) manifests non-local modifications to Einstein's 
equations as quantum corrected equations of motion. That is, at the loop level, matter sources cur- 
vature through a gravitational coupling that is a non-local function of the covariant d'Alembertian. 
If the functional form of the resulting Newton's 'constant' is such that it annihilates very long 
wavelength sources, but reduces to 1/Mp t (M p i being the 4d Planck mass) for all sources with cos- 
mologically observable wavelengths, we would have a complimentary realization of the degravitation 
paradigm- a realization through which its non-linear completion and the corresponding modified 
Bianchi identities are readily understood. We proceed to consider various theories whose coupling 
to gravity may a priori induce non-trivial renormalizations of Newton's constant in the IR, and 
arrive at a class of non-local effective actions which yield a suitably degravitating filter function for 
Newton's constant upon subsequently being integrated out. We motivate this class of non-local the- 
ories through several considerations, discuss open issues, future directions, the inevitable question 
of scheme dependence in semi-classical gravitational calculations and comment on connections with 
■ other meditations in the literature on relaxing of the cosmological constant semi-classically. 

(N ■ I. PRELIMINARY REMARKS 

> 

o . 

y—{ t As is widely appreciated, the cosmological constant (CC) problem [1| is actually three devilishly difficult problems 
rolled into one. To wit, one has to explain why its natural value is not of order unity in Planck units (the so called bare 
cosmological constant problem) as well as explaining why it is not zero (to one part in 10~ 120 ) as current observations 
, suggest [2f] j in addition to explaining why it has only started to dominate the energy budget of the universe recently 
' (the coincidence problem). 

Perhaps the simplest manner in which to appreciate the logical independence of the three cosmological constant 
problems is think of them in terms of their consequences for the lifetime and the structure of the universe. The 
bare cosmological constant problem, posed differently, is simply asking the question why the universe evolves at all 
rather than staying deep within the quantum gravity regime, composed of causally disconnected regions with horizons 
comparable to the Planck size. To ask why it is not zero, is to ask why the universe asymptotes to a de Sitter phase 
in the future. The coincidence problem is to ask why the observed transition from matter dominated evolution to the 
de Sitter phase happens when it does, so that several cycles of stellar creation and destruction have taken place (a 
process that takes of the order of billions of years), so that the heavier elements and the simple molecules that are a 
pre- requisite for the evolution of complex organisms have been synthesized (42j . 

At present, two broad categories of approaches towards resolving the cosmological constant problem can be distin- 
guished. On the one hand there is the anthropic approach (see Q for a partial survey), and on the other, a varied set 
of musings in the literature on the possible dynamical resolution of the cosmological constant problem (see Q @ for 
a review). In particular, it seems reasonable to assume that any attempt that does not resort to anthropic reasoning 
would have be dynamical in nature, especially in considering the coincidence problem which appears to demand some 
form of tracking with matter or some form of relaxation mechanism to alleviate it. We do not wish to delve into 
the relative merits of either approach in this report, whose focus is simply to present a candidate solution to the 
CC problems that arises in semi-classical gravity. This approach- semi-classical degravitation- also qualifies as a 
dynamical approach to solving the CC problems, and the realization that we wish to present in this paper will turn 
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out to have overlaps with several other approaches [4j|. We note in passing that in principle, it would be possible 
to immediately distinguish between any putative anthropic or dynamical resolution to the CC problems through any 
signatures of a time dependent component to dark energy. Such an observation would strongly disfavour anthropic 
scenarios (see Q for interesting indications in this regard). 

We commence this report by reviewing a mechanism which has been known for some time, but whose exact 
realization and exact interpretation when taken as a non-local modification to Einstein's equations remain open 
questions- degravitation [|| |9| , when taken as a phenomenological paradigm, attempts to solve the bare cosmological 
constant (CC) problem directly by nullifying the gravitational effects of the bare CC no matter how large its value in 
Planck units. In doing so it also contains mechanisms to explain what we observe as dark energy in the present epoch 
as a memory effect, in addition to considerably alleviating the coincidence problem [To| . The manner it achieves all of 
this is to posit that Newton's constant is not so much a universal constant so much as it is a scale dependent coupling 
such that its strength for a source with a characteristic wavelength A satisfies: 

MpiG]y(\, L) — > A»L, (1) 
M%G N (\,L)-> 1 A<£. 

where L is some IR length scale (typically taken to be much larger than any cosmologically accessible scale) , and M p i 
is the usual four dimensional Planck mass. Specifically, Newton's constant is to be thought of as a filter function 

G N -► G N {L 2 U) (2) 

that satisfies ([1]), and acts on sources of energy- momentum. In it was shown that linear perturbations away from 
Minkowski space for massive, or resonantly massive gravitational theories coupled to conserved sources results in the 
parametrized filter function: 

8itG n -> , 0<a<l, (3) 

where a = parametrizes massive gravity and a = 1/2 corresponds to the corrections to the brane graviton propagator 
in a DGP-like setup. The bounds on a in the above follow from unitarity [12]. In the case of massive gravity perturbed 
around Minkowski space, such a filter function arises when we integrate out the extra polarizations of the graviton, 
such that the spin 2 part of the metric tensor couples to conserved sources with the resulting filter function (EH) (to be 
understood as an expression only perturbatively). This work was subsequently followed up and generalized in [l3l|[l4j 
to address some of the open issues of the constructions considered in [Til ] . Further developments seem to imply natural 
realizations of the degravitation paradigm (though with drastically different functional forms for the resulting filter 



functions) in higher co-dimensional embeddings 15] If 

We also observe that, although designed to nullify the bare CC by default, the non-local nature of any suitably 
degravitating filter function ([3])- provided we can take its non-linear completion and its existence around a wide 
class of backgrounds seriously- also contains mechanisms to explain the magnitude of extremely small value of the 
apparent CC today as well as considerably alleviating the coincidence problem [lo| : a review of these aspects of 
degravitation is offered in appendix A. As such, degravitation as a paradigm offers a candidate mechanism which has 
the potential to address aspects of each of the cosmological constant problems through only the simple assumption 
that the gravitational coupling Gn vanishes for asymptotically long wavelength sources. 

Motivated by this, we investigate in this report the possibility that degravitation might be nothing other than the 
configuration space manifestation of the running of gravitational couplings (see (l7j for an earlier application of this 
observation). In this way, its existence around arbitrary backgrounds, its non-linear completion and the modified 
Bianchi identities implied by are readily understood, and are precisely analogous to well understood phenomenon 
in conventional gauge theories. That is eq. ([3]), and generalizations thereof, are to be thought of as arising from 
renormalization group (RG) flow. 

We begin this report by reviewing how running couplings can imply non-local equations of motion as a means 
to encode semi-classical physics in gauge theories. We review how familiar quantum phenomenon such as vacuum 
polarization and loop corrections to effective potentials can be viewed as nothing but non-local modifications to the 
equations of motion satisfied by the sources and fields at hand. We consider the same situation when we bring 
gravity in to the picture, later motivating ourselves with a suggestive 2-dimensional example. We also discuss how the 
modified Bianchi identities arc to be understood in this context. Furthermore, in considering the modifications to the 
graviton propagator implied by renormalizations of the form ([3]) , we understand through the spectral representation 
of the renormalized propagator that degravitation, thus realized, is nothing other than the result of having either 
generated a mass gap for gravity under RG flow (as would be the case for a filter function of the massive gravity type 
(a = 0)), or the generation of an infinite tower of massive resonances [Hj][44j]. Having understood the contours of the 
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program we are interested in executing, we then outline our calculational approach to the general problem, paying 
attention to its potential pitfalls and open issues. 

We then proceed to review how gravitational couplings are renormalized upon integrating out various matter 
contents using standard heat kernel techniques. Realizing that any field content that induces non-trivial running of 
Gn deep into the IR would itself have to manifest non-trivial physics in the IR, we then generalize the heat kernel 
technique to incorporate field theories with non-standard kinetic terms and consider various classes of theories that 
feature novel IR physics. We first orient ourselves with examples corresponding to DGP inspired setups and non- 
commutative field theories, among others. Drawing lessons from these examples we are led to a class of non-local 
field theories that formally exhibit degravitation and nullifies the cosmological constant. We motivate this class of 
non-local field theories through several contexts in which it could arise. In a particular scheme, we infer degravitation 
of the massive gravity type (with a = in ©), but where now the non-linear completion of the theory and the 
implied modified Bianchi identities are readily understood. We temper this result with a discussion of the inevitable 
scheme dependence that exists in any such calculation. Although the existence of a degravitating solution appears 
to be a scheme independent conclusion, the precise functional form of the resulting filter function appears not to be. 
We view this as an inevitable consequence of the limitations of our calculational techniques rather than an intrinsic 
aspect of the general program we outline in this report- that non-trivial RG flow in the IR results in a non-local filter 
function as the appropriate gravitational coupling ([2]) . Therefore the perspective of the results contained in this paper 
is that we present the basis of a novel realization of the degravitation paradigm (through which the cosmological 
constant problems are naturally addressed), through which its non- linear completion and existence around arbitrary 
backgrounds are readily understood, and furnish preliminary evidence for this claim. 

II. FILTERING GRAVITY BY RUNNING G N 

In a typical gauge theory, it has been observed that the result of loop corrections in position space, is to effect the 
replacement for the coupling constants [l9j : 

a -> c*o(m) + c ln[-D/n 2 }, (4) 

where /i is some (experimentally determined) reference energy scale (such that a([i) = ao), and the coefficient c 
depends on the renormalization group equations of the theory at hand. In fact, it is relatively straightforward to infer 
from the basic properties of Fourier transforms, that since the bare coupling constant multiplies a product of local 
operators in the bare action, the effects of RG running would be to induce some form of non-locality in the vertex 
at the level of the semi-classical spacetime action once we transform back from momentum space [20(. That is, the 
net effect of quantum corrections to a physical system can be encoded by non-local modifications to the low energy 
effective action, which is then used to extract the correct semi-classical physics. This situation is certainly familiar to 
the reader, as we illustrate with an example from QED, where the interaction Lagrangian reads: 

Ant = eV-fA^x). (5) 

In Fourier space, it is clear that if our coupling runs as a function of the scale (set by the external momentum p) of 
the process at hand (i.e. e — > e(p 2 / /_t 2 )) , then we can just as easily rewrite the coupling as e(D//j 2 ) (where the delta 
function constraint on the vertex results in □ — > p 2 ). This suggests that in position space, we can encode quantum 
effects at the level of the action through the substitution: 

A„t = e(U/p 2 )^A^{x). (6) 

To flesh this out, consider that in QED, it is straightforward to show that corrections to the photon propagator from 
virtual fermion loops results in the running pl| 

a^a(k 2 ) = a \\ + £-ln ^ + 0(e^)l ;k 2 /ml»l, (7) 

L o7T 771 g J 

where a is the low energy fine structure constant a rs I/I37. This implies the non-local form factor (UJ in configuration 
space: 



a 



l + ^-l n [-n/m 2 } + 0(et) , (M 
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Indeed, we know that something of the type has to in fact happen at the semi-classical level if we are to obtain the 
correct corrections to the classical equations of motion from quantum effects. In QED, we know that the electrostatic 
potential felt by an electron due to the presence of an infinitely heavy point charge (with momentum transfer k^ = 
(0, k)) is given by the Fourier transform of the scattering amplitude: 

If we took the fine structure were truly constant (4ira = eg), then the above evaluates to the usual classical answer: 

V(r) = (10) 
However once we account for the fact that couplings do in fact run, we realize that there must be a deviation from 



Coulomb's law once quantum effects are incorporated. It can be shown [211 ], that such a deviation can be accounted 



for through the fact that quantum mechanically, the electron effectively gets smeared out into the charge distribution 

/r1 3 k 
—e^4na(-P), (11) 

which is accounted for in (J9)) through the prescription ((4]): 



f H 3 k 1 

V (r) = j^a(-0)e^ ¥ , (12) 

so that it is straightforward to compute the net effect of the running of a (JJJ to be to give us the quantum corrected 
potential in the limit of a vanishing electron mass (i.e. for mr C 1): 



^ ^ 4ttt 

where 7 is the Euler constant (—27 s» —1.154). Through this example, we specifically see how quantum effects alter 
the classical equations of motion, such that the original QED interaction Lagrangian has to be altered in the manner 
suggested by ^ to capture semi-classically, the effects of vacuum polarization. Specifically, one could either view this 
as the Gauss law (which in the classical equations of motion, results in the potential (|10[) for a point charge) being 
modified to reflect (1131) as the correct potential experienced by an electron, or equivalently the smearing out of the 
electron into a charge distribution given by (fTTj) . Either perspective follows from In a similar fashion, we are 
interested in computing semi-classical corrections to the gravitational action, such that the degravitation paradigm 
might find a realization as a semi-classical effect with ^ being the resulting transcription of In doing so, we find 
that we find a natural accounting for the modified Bianchi identities implied by promoting Newton's constant to a 
filter function- an issue which presently deserves a quick digression. 
Considering the modified Einstein equations implied by © [9] 

Gi = ^G N {L 2 U)Tj}, (14) 

one immediately encounters the issue, that since in general [V M , □] 7^ 0, either the Bianchi identities are not satisfied 
in the above, or new constraints are imposed on a covariantly conserved energy momentum tensor (recalling that any 
covariantly constructed field theory will yield a covariantly conserved energy momentum tensor) . Although this might 
cause some readers anxiety, our perspective here is that the correct Bianchi identity is in fact V AI [G(L 2 D)T^ 1 ] = 0, 
and that this is to be expected considering the fact that in reality, the equations (|14p are actually a paraphrasing of 
the semi-classical equations of motion: 

Gt=^2^). (15) 
m P i 

That is, one obtains the energy momentum tensor by taking the variational derivative with respect to the metric 
tensor of the effective action (rather than the classical action): 



a 1 
— I In- 

^7T V 7 
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(13) 



(16) 
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Where in the example of having integrated out a scalar degree of freedom |2J 



W = / <fixy/-g{x) C eff (x) (17) 



i 



oo 

C eJ f = ^ lim / dm 2 G% s (x,x';fh 2 ), (18) 

with G nc(i , x' '; m 2 ) being the Schwinger-De Witt representation of the Feynman propagator of a scalar field with 
mass m [23| . As we have seen in our previous example, just as loop corrections induce non-local operators (j4} in place 
of coupling constants, which effects the intuitive idea that quantum effects 'smear' out classical sources of fields (fTTj) 
or provide quantum corrections to the classical equations of motion, so the running of Newton's constant will result 
in (fl4l) as a paraphrasing of (fT5)) (45j . Therefore, it should be clear that the quantum corrected equations of motion 
derived from the effective action (£ e ff) are implied by V M [G(L 2 D)T^] = 0, as V^T^ 1 = implies only the classical 
equations of motion, derived from the classical Lagrangian density, L c . 



A. An example in 2 + e dimensions 

In pursuit of further motivation, we discuss a toy example [22| which arises in the context of 2-d gravity- where 
gravity is a renormalizable theory. Duff and Christensen |24j | have computed the one loop beta function when gravity 
is coupled to n s massless spin s particles, and N s massive spin s particles in D = 2 + e dimensions as: 



/?(G) = (d-2)G-/?oG 2 -..., (19) 

where /3 is given by 

A) = |[1 - n 3/2 + ni/2 - no - JVi + N 1/2 - N ], (20) 

From ([TO]) we can immediately read off the existence of a UV fixed point for gravity at G c = e/flo, if /?o > 0. The 
differential equation defining the beta function: 



is readily solved to yield the scale dependence of Gat as: 

a 

1± (m 2 /fc 2 ) 



G ( fc2 ) = ; ■ (») 



where m 2 is constant of integration, determined at the energy scale where we experimentally measurre Gn through 
some local experiment. In configuration space, this would result in the promotion of Newton's constant to the filter 
function: 

G (°) = i, r G 2 /n ^ ' (23) 
1 ± (m J /U) 2 

which suggestively compares with the family of filter functions proposed in [ll| (j3]): 

g ( d ) - i +( J;d)i-. - (24) 

In what follows, we set about the endeavour of obtaining an appropriately degravitating filter function for gravity in 
4-dimensions. Before we do so, we digress momentarily on the physical interpretation of a running Newton's constant. 



B. Generating mass gaps and poles under RG flow 



As discussed in 12j and emphasized in the context of degravitation in 11], any theory that induces an effective 
momentum dependence for Gn necessarily looks like a theory that propagates a massive graviton, or many massive 
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resonances. This is most easily seen around a Minkowski background in the context of a hlter function of the form 
([5]) , which would result in the scalar part of the graviton propagator assuming the form 

A ^ = w+hwr (25) 

where the precise functional form of M 2 (fc 2 ) depends on the value a. We can immediately perform a spectral 
decomposition of the propagator 

A <* 2 > = f *" 2 ra <*> 

where, for example, DGP-like propagators (with a — 1/2 in ((3])) admits the spectral representation [ll| 

2M f°° dm 2 1 , 2 , 2M , . 

P( m ) = 7TH , k > ( 27 ) 



7T 



i 2 + M 2 m 2 -D' PV ; ir(M 2 +m 2 )'- 



and the propagator for massive gravity admits the simple spectral representation: 



1 f°° 7 S(m 2 ~M 2 ) , , 

' dm 2 1 - ■ > . (28) 



-□ + M 2 J -□ 

In the case corresponding to a DGP propagator, we see that the spectral decomposition implies that the theory 
propagates an continuum of massive resonances. For the case corresponding to massive gravity, there is only one 
isolated pole at M 2 , as expected. Evidently any momentum dependence of Gn that admits a spectral representation 
(|2T>1) forces us to interpret the resulting theory as propagating massive resonances. 

Therefore if the effect of RG flow in the IR is to induce an effective propagator of the form (|2"5)l generalized to 
arbitrary backgrounds, evidently one of two phenomena are occurring. Either the pole at zero momenta remains 
fixed and extra poles have been generated by renormalization, or in the case of a resulting filter function of the 
massive gravity type, the pole at p 2 — has been lifted and an effective mass gap has been generated in the graviton 
spectrum. In the former case, any extra isolated poles would have to correspond to new spin states unless there 
were an infinite tower of them. This is due to the fact that any propagator with a finite number of spin two poles 
would correspond to a theory with a finite number of higher order curvature invariants which generically propagate 
ghosts and correspond to unphysical theories. Any physical trajectories in theory space induced by RG flow would 
have to flow along trajectories that do not admit propagating low energy ghosts. In particular, were they to admit 
a parametrization for the renormalized Newton's constant that corresponds to ([3]), they would have to satisfy the 
unitarity bounds < a < 1 [46| . However, the gain in generating ([3]) through renormalization rather than from any 
background specific construction, is that its non-linear completion and consistency is guaranteed provided we begin 
our flow from a consistent (i.e. ghost free) point in theory space. Thus the approach that we espouse in this report 
considers Einstein Hilbert gravity coupled to various matter fields (that describe consistent theories which manifest 
novel IR physics), and attempts to study the resulting RG flow of the coupled system in the spirit of pursuing a 
proof by construction. In this way, degravitation would find a realization as an effect in semi-classical gravity. One 
particular example that we find results in degravitation is obtained by inferring the RG evolution of gravity coupled 
to a toy scalar field with the non-local action: 



S = 



(29) 



which in spite of its unusual appearances, arises naturally in several contexts and defines consistent propagation of 
massless quanta with non-trivial boundary conditions as we discuss further. Before we commence with the thrust of 
our investigation, we take a moment to outline the strategy we will pursue. 



III. THE OUTLINE 



The basic question we are attempting to answer in this report, is whether or not gravitational couplings can be 
engineered to run deep in the IR such that the combination GjvA flows to zero, and that Newton's constant assumes 
the form ([2]) in its action on all other sources of energy momentum. The discussion until now has demonstrated the 
utility and the desirable phenomenology of such a possibility, and the remainder of this report focusscs on attempting 
to realize this possibility in a proof by construction. 
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Clearly we are interested in computing the beta functions for the gravitational couplings in the Einstein Hilbert 
action, and inferring the renormalization group flow that results. There are several independent ways one can attempt 
this calculation, and in this report we focus on the particular approach offered by heat kernel techniques. This is not 
the only way to proceed, and any conclusions drawn in this investigation would ultimately have to be cross checked 
with other available methods. The latter point assumes further importance when we realize that computing the 
running of couplings in a non-renormalizable theory (with dimensionful couplings at the lowest order) is bound to 
involve an irreducible degree of scheme dependence. Although the fact that degravitation is indeed effected in the 
example we consider will turn out to be independent of the scheme, the precise functional form of the resulting filter 
function inherits this scheme dependence. 

Therefore we approach this study in the spirit of delineating the idea that degravitation has a rather natural 
realization as a loop effect in semi-classical gravity, and with a particular technique, coupled to a particular field 
content, furnish evidence to this effect. Any firm claims in the affirmative, would certainly have to be examined with 
other techniques (lattice simulations, functional renormalization group techniques etc). With these caveats in mind, 
we begin our investigation with a review of the heat kernel renormalization technique in curved spacetimes. 



IV. EFFECTIVE ACTION FOR GRAVITY COUPLED TO A SCALAR FIELD 



In this section, we review the derivation of the effective action for gravity when coupled to a minimally coupled 
scalar field employing heat kernel techniques. Starting with the action: 



S 



d A x 



-9 



one can subsequently compute the following loop corrections to the graviton propagator: 

iG R iG B 



k 2 



iG B r.^JG B 



k 2 



k 2 



(30) 



(31) 



where Gb is the bare coupling multiplying the Einstein Hilbert action, and E is given by evaluating the appropriate 
loop integrals, and (ignoring tensor structure) is of the form 



E = 



167T 



;k 2 A 2 , 



(32) 



where c is a number of order unity, and A is the cut-off in our loop integrals. This allows us to infer the renormalized 
Newton's constant as: 



1 

Gr 



1 

Gb 



16tt 2 



rA 2 



(33) 



The net effect of explicitly integrating out canonical matter fields using heat kernel techniques [25| |26J as we shall see 
shortly, is to yield the Wilsonian running: 



G{n) G(/ire/) 12tt 



(n +7ii/2 -4m), 



(34) 



where n is the number of real massless scalars, m/2i the number of Wcyl fcrmions, n 1; the number of gauge bosons 
and /i the running energy scale. One infers from this, the running Newton's constant: 



G( M 2 ) 



l + cG(tf ef )v 2 ' 



(35) 



Although we see immediately that this form is of little use to us in terms of degravitation, the implications of the 
above have been explored in the past to see whether or not the scale of quantum gravity can be brought down to the 



threshold of current accelerator energies by invoking many solely gravitationally coupled species 27] [28] 29] . To recap 



the essentials of the derivation of (l3"5"j) from functional methods, we first consider the effective action for a scalar field 
minimally coupled to gravity: 



(36) 
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Defining the heat kernel: 

H(t) = Tr e- rA = e ~ TXt > ( 37 ) 

i 

with Aj defined as the eigenvalues of A = — □ + to 2 . The effective action can then be obtained via the functional 
identity: 



\Y. lnX ^-\[ 2 dT ^ ( 38 ) 



i-i -i poo JJ(t\ 

W = -lndetA= - 
2 2 



where we are required to regulate the divergent integral with the cut-off in the lower bound of the integral, with the 
definition of the heat kernel as: 

H(t) = [ d 4 x G(x,x;r) = f d 4 x (x| e - T( - n+m2) 0(r)|x), (39) 



where the theta function is necessary as the operator in the exponent has eigenvalues that are unbounded from above 
in Euclidean space. We thus see that G(x,x';t) satisfies the differential equation 

(d T -D x + m 2 )G(x, x'; t) = ^5{t)5\x, x'), (40) 

which justifies the nomenclature. In flat space, for a massless scalar field, one has the usual result 

G{x,x' ] T) = —^e- lB ^. (41) 

In curved spacetimes, we have to expand in terms of curvatures to obtain the result (see e.g. [26j][4 



1 



T 



G 



ff M = i^V2 I dV+- I R + 0(t\R<) . (42) 



A derivation of this result in a different context will be presented shortly. For our immediate purposes, we note from 
the above and from (f38|) . that the renormalized Newton's constant is given by: 



1 1 1 

Gr G re f 12ire z 

One can infer the Wilsonian running of Newton's constant by imposing an IR cut-off [48[ : 

w= .if' iT m, (44) 

2 Je 2 T 

and computing the renormalized coupling at a reference energy scale, relative to which we compute the renormalized 
couplings at any other scale, which we can then write in a manner independent of e as: 

1 1 a 2 

- C 7^Z> (45) 



G(/z 2 ) G ref 12n 

where all the e dependence has been swept into the (renormalized) quantity G re / which we determine experimentally, 
and where c is some numerical factor that depends on the number of scalar, fermionic and gauge field degrees of 
freedom. Similarly, one infers the renormalized cosmological constant: 

A R = A B + (N b -N f )^. (46) 

If our goal is to engineer a running for Newton's constant such that the gravitational coupling it implies vanishes in 
the far IR, we encounter the immediate realization that conventional matter and gauge fields do not do the job. Given 
that any direct measurement of Gat (say by Cavendish type experiments) sets its scale to be l/M 2 t at laboratory 
scales, we require that any running induced by matter or gauge field couplings be such that it not only has the right 
sign, but is also strong enough to render gravity a dud force on extremely large scales. One can easily see through (|34[) 
and (|35l) that no massless scalar, fermionic or gauge degree of freedom can induce a suitably degravitating running of 
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Newton's constant. The same can easily be shown to be true of massive degrees of freedom as well. These observations 
are perhaps painfully obvious, as why should renormalization of these fields- a process which involves the subtraction 
of infinities due to divergent loop integrals in these theories- run deep into the IR? The answer is that it does not, 
unless the theory itself presents us novel physics (RG flow) in the IR. This novel physics would have to be such that 
were we to somehow obtain the beta function for this theory coupled to gravity, it would integrate to cause M^Gjv 
to run from in the far infra red, to 1 at all cosmologically accessible scales. 

Phrased in this manner, we are simply asking the question of whether or not it is possible to engineer a field content 
such that gravity has a positive beta function from the far infra red to Hubble scales, and such that M^Gn —> and 
GatA — > in the far infra red. Although the answer to this question in all generality is beyond our means at the 
moment, we attempt this investigation in the spirit of pursuing a proof of concept by construction. In attempting 
to understand the class of theories that degravitate by construction, we also uncover hints of the possible underlying 
physical mechanisms for this non-trivial RG trajectory in the IR which we will elaborate on further. Before we do so, 
we observe the general features of a theory that would serve our purposes. 



V. RUNNING G N IN THE IR 



In curved spacetimes as we shall see shortly, one can compute the curvature expansion of the heat kernel to find 
the general form: 

H(t) = I d A x G{x, x; t) = A(r) [ dV + g(r) [ dV R + Q{R 2 ). (47) 



In order to obtain the effective action, one integrates the above as: 

W-\r*®., (48, 

where the cutoff in the lower limit in the integral (j48]) is necessary to regulate the integral. Noting that the (bare) 
Einstein-Hilbert action is of the form: 

' f dV R+ -5- / dV A B , (49) 



16nG B J 8tt 
we see that the renormalized gravitational couplings are given by 



■k—k- + »*r 9j r> (50) 



A/, = A 1; - 4tt / (51) 

'l T 



To obtain the Wilsonian running of Newton's constant, we compute the expression: 



1 1 ..ftrf ^, (52) 



G(p) G B 

and expressing Gb in terms of some reference G(/i r ) (and hence absorbing the usual infinite renormalizations required 
by UV divergent loop integrals) at some reference (UV) scale [i r via 

1 - 1 8./""^, (53) 



Gb G((i r ) J e 2 t 
we can substitute back into (|52[) to compute the running of Newton's constant as: 

Gjv(/i) = GnM 3 „ . (54) 

The scale [i r is necessarily a UV scale (/x~ 2 < /i~ 2 ), as we can only directly measure the strength of Gjv via Cavendish 
type of experiments, and as laboratory scale experiments, are in the far UV as far as cosmological and astrophysical 
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scales are concerned [49(. We also note here, that in order for degravitation to have meaningfully achieved its 
main purpose, in actuality we require that the product G^(/i)Afl(/i) vanish in the IR. That is, given that whatever 
makes Newton's constant run will also make the cosmological constant run in the IR, in order to degravitate the 
bare cosmological constant, we require that the product GjyA (which is what would source curvature) vanish in 
appropriately dimensionless units in the IR. In this way, the goals of the endeavour of this paper might not seem 
to be so different from the many approaches to relaxing the cosmological constant via RG running that exist in the 
literature [301 ] . However, as we have noted in the introduction and in the appendix, solving the bare CC problem 
(by construction) is but one aspect of degravitation, which is a more general mechanism. Considering G» as a 
filter function acting on all other sources of energy momentum (including CC like contributions), has advantageous 
features for addressing the other two cosmological constant problems, as well as cosmology in general [l(| due to the 
intrinsically dynamical nature aspect of this prescription. 

The perspective one should have then of the sequel is as follows: we consider the effects of integrating out a particular 
scalar field in some hidden (dark) sector (i.e. this field couples only minimally couples to gravity and nothing else). 
We are then in fact, obliged to integrate out this field when considering the low energy interactions of all other fields 
coupled to gravity. Upon integrating out this field, we will have satisfied our objectives if GatA/M 2 ; runs to zero, and 
the resulting Gn becomes a filter function when acting on all other sources of energy density (treated classically). 
In this way, any other contributions to the cosmological constant are also dynamically rendered gravitationally inert 
in such a manner that can also address the other two cosmological constant problems (see appendix A). We orient 
ourselves first with two familiar examples. 



VI. SCALAR THEORIES WITH IR MODIFIED PROPAGATORS 



In considering which theories when coupled to gravity might render non-trivial RG flows for the gravitational 
couplings in the IR, it is a fair assumption to draw upon theories which in and of themselves exhibit, or reflect novel 
physics in the IR. One such example is furnished by a minimally coupled a scalar field with the Euclidean signature 
propagator: 



1 



□ 



(55) 



Such a modified propagator naturally arises on a DGP brane (3lj as a result of accounting for bulk propagation in 
the effective theory (with an explicit brane mass term included). We generalize the standard heat kernel technique 
for computing renormalized couplings induced by such a modified propagator in appendix B. There, we derive that 
the function A(t) which determines the renormalization of the cosmological constant through (|47[) is given by (|125p : 



A(r) 



16tt 2 



1 M 

^ + 17™- 



ttt(6 + M 2 t)[1 =f Erf(My/r/2)] 



8t 2 



(56) 



where the signs ± correspond to the choice of the propagator (O^pMy/— □ — m 2 ) -1 , with the lower sign corresponding 
to the so-called 'self- accelerated branch'. Correspondingly, g(r), which determines the renormalization of Gn through 
(jnnj) is given by: 



9(r) 



16tt 2 



1 M 2 
7^ + -F^T 



M 



(l t £r/[MV7/2]) (3 



M 2 



96 1 32t 

We note from (f54|) . that in general, in order to achieve the running G(/x)/G(/i r ) 
/i — > 0), we require that the integrand <?(t)/t decay no faster than 1/r at large r 



(i 



G(/z)M 2 - 
From ([57]). 

renormalization of Gn at IR scales vanishes at large r for the choice of the propagator (□ — My 
inferred from the asymptotic formula for large r [321 ]: 



1 - Er/[Mvr/2] 



,-M 2 t/4 



irtM 



1 



Q{t- 



(57) 

-> in the IR (as 
we note that the 
— D) -1 , as can be 



(58) 



The so-called self-accelerated branch however, formally exhibits the requisite behaviour in the IR. When we tune the 
brane mass m of our model to be such that m 2 = M 2 /4 (the choice of which we elaborate upon shortly), the dominant 
IR (large r) behaviour of ([57| goes as M 3 /96 • y/nr/Wn 2 . This would result, after obtaining the effective action via 
and (|49[) in a filter function for long wavelength sources of the form: 



G N {M- 2 D) = 



Gn(Ht) 



1 + y/irG N (fi r ) 



M 2 



M 



M 

fly 



(59) 
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By defining G as 



we can rewrite the above as 



G N (M- 2 D) 



Gn{^t) 



1 - 



Gn(/Jt) 



M 



M 



(60) 



(61) 



We note with interest how the above compares to © for the case a = 1/2 [Tlj . We can also compute how the 
cosmological constant runs via (|125l) as: 



with 



A(m) = A(/Xr) 



A(// r ) = A(/i r ) 



M 3 _ 



m 3 _ 



(62) 



(63) 



recalling that the infinite renormalizations of the CC coming from the standard UV divergences are all contained in 
A(ju r ). From ((61]) and ([62]l. it is clear that G(^)A(^)/M 2 , -> in the IR, with the degravitating filter function (|6Tj) 
resulting upon all other sources. 

In spite of its appearances however, this example is unphysical. We observe in this that the physical origins of this 
apparent degravitation lie firmly in a new class of IR divergent loop diagrams. Upon cavalierly computing the one loop 
effective action given the form of the DGP propagator on the self-accelerated branch, we find the loops responsible for 
the Wilsonian running that appears to cause the vanishing of Newton's constant in the IR are precisely those diagrams 
that correspond to poles on the positive real axis in momentum space. In particular, consider one of the diagrams 
that contributes to the one loop effective action with external graviton legs and internal scalar field momentum k: 



AT, 



d 4 k 



1 



k 2 ±Mk- 



(64) 



It is clear that we can write the denominator as (k — 0Ji){k — u>%), with the u>i being real and positive roots on the self 
accelerated branch. The choice m 2 = M 2 /4 makes the denominator take the form of a perfect square (k — M/2) 2 , 
preventing us of regularizing the integral via any sort of principal value prescription. Therefore, we understand this 
apparent degravitation as arising from IR divergent loops due to states with the dispersion relation k 2 = M 2 /4. 
Although these states correspond to plane wave solutions of our field equations in Euclidean spacetime, they are 
exponentially unstable when continued back to Lorentzian signature and do not correspond to a sensible theory, i.e. 
they describe tachyons. The self-accelerated branch is unstable precisely because of these tachyonic modes [33| and 
therefore does not admit a meaningful Minkowski space interpretation. We therefore discard this example, but glean 
from it the diagrammatic roots underlying the running of Gat in the IR for this type of propagator. 

We now turn to another example which manifests novel physics in the IR to gain further traction on our problem- 
a class of widely studied (non-local) field theories whose one loop renormalization of gravitational couplings can also 
be computed. We consider the propagator: 



1 



□ - 



1 2 ■ 

— m z 



(65) 



FLJ 



where 6 has dimensions of length 2 . This propagator can be thought of as arising from loop corrections to a non- 
commutative scalar field theory in Euclidean spacetime [34|[35|. or it can be thought of as the free field limit of a 
translationally invariant, renormalizable non-commutative field theory [36j]. As before, we are not interested so much 
in the specific origins of this theory in so much as shall use it as an theoretical laboratory in which test the ideas 
which we are interested in [Ho} . It is straightforward to calculate the renormalizations induced by this example (using 
(jl39[) in appendix B, and generalizing to the propagator (|65jl ). We find straightforwardly that 



A(t) = I d x G(x,x;t) 



8tt 2 



-K 2 [2t/6] 



(66) 
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with A(r) again defined as in (l5Tj) renormalizing the cosmological constant, and with K2[x] the modified Bessel 
function of second order. Similarly, Newton's constant is renormalized as 

air) = 3 [^[2r/0] T0 K 3 [2r/6] + ^K^r/O}] (67) 



-K 2 [2t/6] 



8tt 2 2 6 

with g(r) defined as in (|50|) . and the second line follows from the first via a recursion identity for the modified Bessel 
functions 32]. Once more, we find that this model does not exhibit degravitation. This is because the asymptotic 
form of the modified Bessel functions is exponentially decaying [32| : 

lim K n (z) ~ J^-e-\ (68) 

and therefore g(r) / r decays too fast to cause the vanishing of Gat in the deep IR. However by deliberately introducing 
a tachyonic instability as above, one could again formally induce a degravitation filter, but this time of the massive 
gravity type (eq. ([3]) with a = 1) [5l|. It seems neccesary therefore that in order to consistently effect degravitation, 
one has to consider theories with divergent diagrams that yield the appropriate IR asymptotics for g(r) without 
signalling an instability in the theory at some finite momenta. Moreover this divergence would have to be physical 
in nature rather than an artifact of any particular scheme. We consider a class of theories which posses the desired 
features next. 



VII. DEGRAVITATION- A TOY MODEL 



We reconsider the general expression for g(r), the quantity which determines how Gm is renormalized through ([J 
Beginning with (|1 39[) in appendix B, denoting x — k 2 and performing the angular integrations, we find that 



T 1 f°° T 2 

9( T ) = 7T ~rz — 7t / dx e- Tu{x) xu'(x) + u"{x)'-- - -x 2 u'{xf . (69) 



Realizing that the first two terms in the square brackets sum to (u'x 2 /2)', and upon integrating by parts, realizing 
that the last term is minus one half of the sum of the first two terms, permits us to condense all of the above into the 
expression: 

- iis?r 4 "^"s( ! r) ( ™» 

dx e-^W. 



12 16tt 2 J 



We now first perform the r integral (which corresponds to a change in the renormalization scheme from before), 
resulting in: 

PI) 



12 16tt 2 _ 

where now instead we have to regulate the limits of the x integral (recalling x — k 2 ). Considering the example 
of a canonical scalar field (with u{x) — x) yields the standard UV divergence found in (|3"3"j) . Inserting the two 
different modified propagators considered in the previous section also yields the same UV divergences (~ A 2 which 
were absorbed into G re f), and give us IR divergences at precisely the same energy scales (at the tachyonic poles), 
though with different functional forms around this divergence. Thus the inherent scheme dependence of computations 
in semi-classical gravity rears its head here- a seemingly inescapable state of affairs (see [13] for a discussion). In 
general, in any enterprise involving computing the running of gravitational couplings such as the one undertaken here, 
the only strategy we have to proceed with in ensuring the physical sense of our calculations is to repeat them in various 
different schemes, and at the very least, see if the qualitative nature of the conclusions is scheme independent. We 
presume that the vanishing of A([/,)G(p)/M 2 l in the far infra-red (and not at some finite energy scale) is unambiguous 
enough a feature which one can repeatedly test in different schemes. We bear this important caveat in mind in what 
follows. 
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Resuming from (|71[) . in writing our kinetic term as u(x) = exp[g(x)] with x := fc 2 , we see that we require the 
integrand 

dx [xg'{x)] 2 (72) 

to have a divergent lower limit. This would only be possible if xg' diverged as some negative power of x for small x. 
An example of a non-local scalar field theory which accomplishes this is given by the action: 

S= f d 4 xy/=g~ <j)[e m2 :/a D}<j), (73) 



In spite of its unusual appearances, we motivate this theory shortly Starting with (|71|) . upon regularizing the 
lower limit of the integral with the IR cut-off jj? the generalization of ([5B")) becomes: 



1 1 1 



G(fi) G(ji r ) 24tt J 3 u 



(74) 



where again, the UV divergent contributions have been absorbed into Grfi r ) analogously as between (f50f and (j53"j) . 
For u{x) = xe~ m / x , the generalization of (|54p becomes: 

g (m) = q ¥k rr ( 75 ) 

which impiles the filter function 

G(D/m 2 ) = Q ^r) (76) 

1 + MM — 

which corresponds to a filter function of the massive gravity type, and where we defined: 

G(fj, r ) = ° {Pr) 2 2 . (77) 

We discuss the scale m that appears in the above shortly. The renormalization of the cosmological constant A is given 
by: 

If 1 f r m 2 ] 
Afl = Ag + — I xdx ln[u(x)] = As + — I xdx In x . (78) 

47T J 47T J L x - 

The action (|73|) results in the standard cut-off 4 divergence (easily seen from the fact that the kinetic term asymptotes 
to the standard kinetic term for large momenta), and yields the running in the IR: 

2 

— TYi 

A(/x) = A( Mr ) + — ^ (79) 



with the definition 



2 

m 



A(p r ) := A(/i r ) - ^- M 2 . (80) 

Hence the product G(/i) A(/i) — >• as /1 — > 0, and all other cosmological constant like sources (which will not renormalize 
Gjy in the IR) will be acted upon by the filter function (|76[) . Therefore this theory appears to offer us a semi-classical 
realization of degravitation where instead of Einstein gravity, we have its semi-classical incarnation: 

Gi = ttt K s-T* 1 , ( 81 ) 

JW P' 1 + 247TM 2 , TT 
pi 

where we have implicitly taken the scale \i r to be such that we measure as the gravitational coupling at 

laboratory scale. We note that in order not to affect the behaviour of gravity at sub Hubble scales, we require that 
the mass scale to in the above be such that 

to < KT 30 M pi « 10- 3 eV, (82) 



14 



whose value we return to discuss shortly. More pressingly, we wish to draw caution against the precise form of (|81l) . 
In the scheme where we compute the running within the proper time formalism (|70[) (instead of switching back to 
momentum space as we did in ([71])). we find that although degravitation still results, the precise functional form of 
the running differs. Evaluating the running for Gn that results from (|54|) and using (|T0[) with u(x) — xe~ m ' x results 
in 

- — T \ JT, (83) 



M 2 

and although the integral we have to evaluate to compute the running of the cosmological constant is not possible 
analytically, we find that \{t) as defined in (|4"Tj) numerically evaluates to something that is accurately fit by: 

4 

777 C 

where c ~ 0.2 is a numerical co-efficient. This results in the running: 

4 

A( M ) = A(fx r ) - c^-ln[ln(m 2 /fi 2 )], (85) 

from which we see as before that Gjv(/i)M5 — > and Gjv(/Lt)A(/j) — > as /i — > 0, however the IR asymptotics of the 
running is markedly slower in this scheme. Thus, although we have demonstrated by example that rcnormalization 
group flow in the IR can indeed result in a filter function for Newton's couping ([2]) with all the properties required 
of degravitation, more work remains in order to place this claim on firmer footing. In particular, a full functional 
renormalization group H analysis might be amenable to the context we are interested in. 

We now return to motivating the action (|73[) . We first observe (as is perhaps not widely appreciated) that inverse 
powers of the d'Alembertian are guaranteed to arise in the IR due to finite terms in the effective action for a scalar 
field 3!) , where operators of the form 

Ar ~ ly (86) 

are shown to appear in the IR in interacting massless theories (late times (i.e. large r) in the so-called proper time 
formalism used in our heat kernel calculations). Kinetic operators of the form □e m / n can perhaps then be thought 
of as modelling a non-perturbative resummation of all such corrections, which we would then subsequently integrate 
out entirely to result in the runn ing, we require. Indeed it is precisely operators of this form that could encode, for 
example, particle creation effects [39j as well as vacuum polarization in non-trivial backgrounds- the consequences of 
which for dark energy and the cosmological constant problem are discussed in Therefore in spite of its unusual 

appearances, (I73|) might even be a necessary modification to a QFT interacting with gravity at very long wavelengths. 
This is our preferred interpretation, though we wish to mention another motivation for (|73[) which can be inferred 



from the context of non-commutative field theories. In this context, it is known 35| that the low energy effective 
action receives non-local loop corrections of the form: 



d i x<f>(n + Jgj^, (87) 

where 9 has dimensions length squared and parametrizes the degree of spatial graininess [53] . Expanding the action 
(|73p at finite momenta results in the expansion: 

S= / d 4 x 4>(p+ !^ + m 2 + ...)0, (88) 

from which it is tempting to conclude that (|73[) might again, be some form of a non-perturbative resummation of loop 
corrections in a non-commutative field theory [H3]. We also note that the theory described by ([73| when considered 
as the limiting case of a sequence of regularized theories (see footnote under (|73p ). describes a causal theory with a 
well defined Cauchy problem (40j[55|. 

Returning now to the scale m in the toy action (|73j) . which we see independent of the precise form of the filter 
function f (|76j) or (|83|l ) that results, the requirement that gravity remain unaffected for all sub- Hubble scales requires 
777 < 10~ 3 ey. We immediately see the naturalness of this requirement depends on the particular physics that gives to 
([73]) . however in both contexts described above, it is straightforward to see that this amounts to a technically natural 
tuning. In the context of having resummed all finite corrections to the effective action of some light scalar degree of 
freedom (in the dark sector) , we see that this scale m is by definition the end result of having summed an infinite 
tower of quantum corrections. We now offer our concluding thoughts. 
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VIII. CONCLUSIONS 



In this report, we have offered preliminary evidence that degravitation |8][9j can be rather naturally thought of 
as an epiphenomenon of semi-classical gravity emergent in the IR. In realizing degravitation as a manifestation of 
renormalization group flow of gravitational couplings in the IR, its sense as a non-local modification to Einstein 
gravity, its non-linear completion around arbitrary backgrounds and the question of the modified Bianchi identities 
that result are readily understood. That it not only solves the bare cosmological constant problem, but also offers us 
the appropriate dynamics to address the other two cosmological constant problems should serve as strong motivation 
to further study this mechanism. In particular, that this type of long distance modification to gravity might be forced 
upon us through non-perturbative resummations of terms that contribute finite renormalizations to the effective action 
for a canonical scalar field (39l j , is a tempting idea. There is deep conceptual appeal in the idea that the cosmological 
constant problems might be resolved, rather than aggravated by loop effects in our effective theories. 
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X. APPENDIX A 



In this appendix, we offer an overview of the mechanisms through which degravitation addresses the various cos- 
mological constant problems. Some of the details concerning what follows can be found in pJp| . We assume a filter 
function of the form: 

G N (L 2 D) = -^y, (89) 
1 + TT 

where G^ = V-^pj is the usual laboratory scale 4d Planck mass, and M is a mass scale which defines the IR filter 
scale L (M = m -j^ in the examples derived in the previous sections). Clearly given the spectral decomposition of any 
particular filter function (|26p . we can easily generalize the results of this section to other functional forms. Equation 
is to be understood in the context of the modified Einstein equations: 



G£ = 8ttGn(L 2 0)T^ . (90) 
To avoid the complexity of tensor structure, we consider the trace of both sides of the above equation in what follows: 

R(x) = -8ir[G N (L 2 n)T](x) (91) 

Given (j89j) we can recast the above as 

R + M 2 U- l R = -8nG N T, (92) 

or more concretely: 

R(x) = ~8nG N T(x) - M 2 J dV yJ-g(xf)G°(x, x)R(x'), (93) 
where G°(x,x') is the Green's function for the dAlcmbcrtian operator- 

□*GV^) = 4=L (94) 
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for which we chose the appropriate causal boundary conditions (so that no sources in the future lightcone of a given 
spacetime point affects the geometry at that point). With this choice, the support of the integrand in extends 
only over the backward lightcone of the point x. One can repeatedly iterate the process above to obtain: 

R(x) = -8itG n \t(x) - M 2 (G°,T)(x) + M 4 (G°,(G°,T))(x) - Af 6 (G°, (G°, (G°, T))(x)l + (95) 

where: 



{G°,T){x):= / dMy/-g(at)(?{x,x')T{x'), (96) 



and the subscript on the integral ▲ emphasizes that the domain of integration is the backward lightcone of the point 
x. Because of this stipulation, the nested integrals in the above are ordered integrals, with each integration variable 
in the past lightcone of the one proceeding it. We can now proceed to evaluate (|9"5j) for various sources. 

We commence with a bona fide spacetime zero mode (i.e. the bare cosmological constant). In that case T = 4A, 
and the nested integrals each evaluate to a constant and can straightforwardly result in: 

M 2n (G°,(G°,...,(G Q ,T))(x) =4AM 2 "A™, (97) 

with 



A := / d A x' y^g[x T jG°(x,x'), (98) 



where the domain of integration now stretches infinitely to the past and therefore integrates to (an infinite) constant 
(A — > oo) [56j. Therefore in resumming (|95[) we find: 

R = - 8nGN TTM^A=°- (99) 

Thus the bare cosmological constant decouples from geometry. This is perhaps more directly seen from (1891) , where the 
fact that the d'Alembertian is a self adjoint operator allows us to simply replace the appropriate eigenvalue in place of 
the operator argument in operator functions. Since in any given spacetime the d'Alembertian always admits constant 
functions as zero modes (among others), substituting the null eigenvalue in place of □ in also demonstrates the 
desired result. 

We now consider a spatially homogeneous source that corresponds to a step function in time in a given FRW slicing. 
That is, we consider a source of the form: 

T = —AV6{t — to). (100) 

This source is to caricature the condensing of a potential dominated epoch at some finite time in the universe's past. 
Specifically, since we are only considering the trace of the energy-momentum tensor, one can imagine this ansatz 
to model the onset of inflation from some pre-existing quantum gravity phase, where a radiation dominated epoch 
makes a transition to a potential dominated epoch at some time to- We don't tether what follows to any particular 
scenario, and simply consider (|100[) as a prototypical source. When substituted into (195]) the time ordered nature of 
the integrals allows us to recast the ordered iteration to obtain the degravitated source as (loj |: 

T deg (t) = -Wd(t - t ) [l - M 2 Sl(t) + - M 6 ^ + ...] , (101) 

with Q(t) given by 

fi(t) := J G°, (102) 

with the domain of integration consisting of the past lightcone of the observer at t up to the nucleation event at 
to- Therefore at t = to, the integrals all vanish given their vanishing domain, after which they steadily increase and 
asymptote to an infinitely large value. The series (jlOip sums readily to give us: 

T deg (t) = -W9(t - to)e- M2n ^ . (103) 
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Therefore the curvature sourced by a step function like source (llOOp is given by [Icj : 

R(t) = 8irG N 4Ve~ M2 ^ G °. 



(104) 



In its full generality, this equation is a very complicated non-linear integral equation. However we can infer some 
general features from its form. In particular the argument of the exponential grows as: 



R(t) 



,-M 2 (t-t y 



where t — to is the co-ordinate time that has elapsed since the step function has nucleated [57 
easiest seen by considering a test source around a Minkowski background. In that case, G° r 
/ A g ° ~ (* - f o) 2 , resulting in (|105[) . From this, we can infer the characteristic timescale 

t ~ 1/M. 



(105) 

This is perhaps 
l/27T 2 r 2 , so that 



(106) 



Imagining now a toy universe where some energy density condensed out of a radiation like epoch according to (|100l) . 
how long would we have to wait for this energy density to degravitate to the currently observed density that we 
accord to dark energy? That is, given V of any particular magnitude, how long would we have to wait before this 
energy density decayed so that R/Mi ~ 10~ 120 ? Evidently, this would be (for initially GUT scale energy densities: 



10~ 12 ) when the argument of the exponent becomes ~ 248. Therefore, 

M 2 T 2 ~ 248. 



(107) 



Given that in the examples considered in the previous sections, M = m-gf-, and given that requiring gravity to remain 
unaffected up to present Hubble scale Ho implies (as is easily inferred from (|8Tj) or ([83)) ) that 

< 1, 



9 9 

m m 



M 2 H 2 

so that substituting the definition of M in (|107[) results in 



And so, 



H^T 2 > 248. 



T > 16H n 



(108) 



(109) 



(110) 



Therefore, if we were to degravitate an initial GUT scale energy density all the way down to the currently observed 
energy density ascribed to dark energy, we would have to wait a period approximately 16 times longer than the current 
age of the universe. Clearly this particular example does not plausibly account for what dark energy might be. 

However, the non-local nature of the degravitation filter ([59")) results in some very striking features as far as cosmology 
is concerned. In particular, it implies a memory effect of past energy densities [10] that might account for the presently 
observed component that is causing the expansion of our universe to accelerate. To see this, we first note that we can 
recast the action of ([55)) on some source p as: 



(x) = p{x) - M 2 / d 4 x' v / ~g(x')G(x,x')p(x'), 



(111) 



where now G(x, x') is the Greens function for the operator □ + M 2 , and Pdegrav(x) is the effective source obtained by 
action of the filter, with Gn = l/M 2 pl explicitly factored out as in (|M[). We now consider the effect of a test source, 
homogeneous in spatial extent but a delta function in time. We work in the test approximation around a fixed de 
Sitter background with co-ordinate interval: 



ds 2 = 



1 



H 2 V 2 



{drj 2 — dxidx 1 ). 



We are interesting in evaluating: 



Sdegraviv - Vi) = S(r) - r)i) ~ M 2 drj' ' d 3 ' x' ' V - g(x')G(x , x')6(rj' - Tji). 



(112) 



(113) 
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Since our source is spatially ho mog eneous, given the standard expression for the Green's function for a massive field 
on a de Sitter background (see [1(| for details), we can immediately perform the spatial integrals to yield: 

dr}' j d*x'^W)G{x,x')5{r,- ni ) = J d/-^^^^sinh[vln{r 1 M)}8{r 1 - m) (114) 



2H 2 vrn l(-T]i) 3 / 2 -" (-m) 3/2+v 
given that conformal time and cosmological time are related via 



-Ht 



*> = — fT' (115) 



we can rewrite the above as: 



M 2 

Sdegraviv ~ Vi) = #(v ~ Vi)~E 77- e~ ^ {t ~ u) Sinh[Hv{t - U)l (116) 

Hue nti 



with v defined as: 



Thus the action of the degravitation on a delta function source is to yield a decaying step function, whose magnitude 
immediately after the rji is immediately found to be 

Sdegrav(V - Vi) ~ 0(V ~ T)i)M 2 A(t), (118) 

where A(i) has dimensions of time, but is an otherwise decaying function. Thus the delta function source at time ti 
(which can accurately model the trace of the energy momentum tensor of a universe undergoing a burst of inflation 
from some initial radiation dominated epoch followed by reheating- provided of course the duration of inflation is far 
less that M _1 ), under the effects of the degravitation filter, results in an afterglow energy density which would mimic 
a cosmological constant with a magnitude in Planck units naturally set by the filter scale: 

Ml Ml Ml Ml~ W ' (11J) 



where the latter inequalities follows from the definition M = m-p^, and the requirement (1821) that the filter function 

not affect gravity up till horizon scales. In the context where we take the limit H — > in the above (so that the 
background spacetime in which we place this test source tends to Minkowski spacetime) , one can interpret this as the 
generation of an 'afterglow' energy density that subsequently sources a de Sitter phase. 

Thus not only is the present scale of dark energy readily explainable within our framework as a memory effect, whose 
magnitude finds ready, and technically natural explanation as the filter scale m, the dynamics also implies that the 
co-incidence problem is greatly alleviated. The reason for this is that even though vacuum energy densities do not scale 
with the expanding spacetime, they decay over time as a result of degravitation and consequently, the coincindence 
problem (which obtains its severity from the fact that vacuum energy densities do not redshift) is greatly ameliorated 
[l0( ■ Thus we hope to have demonstrated that in making the simple assumption that long wavelength sources decouple 
from sourcing spacetime curvature, all three cosmological constant problems are readily, and naturally addressed. 

XI. APPENDIX B 

In what follows, we generalize the derivation of the heat kernel (|^|) adapted to propagators with modified kinetic 
terms. For concreteness, we begin with the example of a DGP like propagator (|55p but immediately generalize our 
results to more general functions of □. Our treatment borrows from that reviewed in [41| . Working in the expanded 
basis \x, t), where we augment the usual 4-d position eigenkets \x) with the basis for r space |r), we see that in order 
to compute the heat kernel ([57]) . we need to evaluate: 

G(o; ) y;t) = (a; > *|- r ^ - ] = ■ ? |y,0). (120) 
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We realize that we can expand the operator in the denominator in a curvature expansion around a locally fiat 
coordinatization of our manifold as: 



9t - [□ - NlV^a] + m 2 = d t - [Do - My/-n ] + m 2 }A (121) 

M 



where Dq = d^d^ is the usual flat space d'Alembertian, and where 

SD:=[a' M d li a v + V l d ll ]. (122) 

The inverse of the operator grouping A is simply the Green's functions in flat coordinates that we compute presently. 
By (|40[). we see that we first have to find the appropriate Green's function for the diffusion equation implied by the 
modified propagator in flat coordinates. Writing the heat kernel as G(x, x';t) — G(x, x'; t)9(t), one has the following 
boundary value problem to solve: 

G(x,x';0) = 8 A {x-y) (123) 
6(x,x';t) = (a x - MV^D x -m 2 )G(x,x';T). 



By translational invariance, we Fourier decompose G and obtain the following solution which satisfies the above with 
the appropriate boundary value: 

_ 2 

G(x,y;r) = ^— [ d * ke -iHx-y) e -T[k 2 +Mk] (124) 

(27T) 4 J 



dk 



Jj(fcr) c - T \k*+Mk] k 3 



47r 2 J kr 

Where J\ is the first order Bessel function of the first kind. In the absence of the term proportional to M, one 
straightforwardly recovers (|4ip. When y — x, we can easily evaluate (| 124|) to yield 

e- mV r 1 M 2 „ e^^F(6 + M 2 r)[l T gr/(M V?/2)] i 
G(x,x;t) = — — -J „- + — 5-0 , (125) 



16tt 2 



t 2 ' 4t ' "" 8r 



where the different signs correspond to the choices □ ± My/—D for the kinetic term. We note that implicit in (|4Q]) is 
the fact that we are working in units where the diffusion co-efficient is taken to be unity, thus giving r dimensions of 
inverse mass squared. The function A(t) (which determines the renormalization of A through (|51|1 ) is given by (|125[) . 

Proceeding, we realize that in a locally flat (Riemannian normal) co-ordinatization around the point x, one can 
expand the metric as 

9n\(x') = V»\ - ■^^(aOy'V, (126) 
where y^ = x'^ — x^- allowing us to compute the operator coefficients in operator grouping B as: 

^ = ^ViA ( 127 ) 

b» = -W^. (128) 



Performing the expansion: 



G(x,y;t) = ( x ,t\-J-^\y,0) = (x,t\j + jBj\y,0), (129) 



which is equivalent to a curvature expansion, we can evaluate the curved space modified heat kernel as 

G(x,y;t) = (x,t\A- l \y,0) + [ / d 4 z(x, t\ A~ x \z, t)(z, t\BA~ x \y, 0). (130) 
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Denoting the flat space modified Green's function (1124[)(|125[) as 9(x, x';t), we rewrite the above as 



G(x,y;t)=g(x,y;t) + J dr J d 4 zQ{x, z; t - r) [^"(z) — — + b"(z)— ] (l + —=^Q{z, y; r), (131) 

where first term in the above for x = x' is given by (|125|) . We next consider the second term in the above. We begin 
by inserting ([T2T )) ([T25 ]) . to obtain: 



dr / d A zQ(x, z; t — r) 



-R^ v z p z x — — - -R^z x — 

L3 p x dzudz" 3 8z». 



M 



G(z,y,r) 



We now substitute the Fourier expansions (|124|) to obtain: 



(132) 



dr / cf 



(2tt)* 



d A k I d 4 pe~ ik ' {x ~ z) e -{t-T)[k 2 +Mk] e -r[p 2 +Mp] 



IflM - Z P 2 A A A _ 

3 " A c^dz" 3 9zM 



In order to proceed in a more orderly manner, we introduce the function 

u(s) := s + My/a, 

so that we can express the above as: 



f M \ 



-ip-(z-y) 



(133) 
(134) 

(135) 



9(r) 



dr / d z- 



(2tt) 8 



d A k / d A pe -*k-(x-z) e -(t-T)u(k ) e -ru(p ) 



3 ' A dz^dz" 3 dz^ 



u / (p 2 )e~ ip ' ( *~ ! ' ) ' 



(136) 



where the prime refers to the derivative of (|135|) with respect to its argument, and we identify g(r) as the co-efficient of 
the spacetime integral of the Ricci scalar as defined in (|47[) . In this form, we can straightforwardly consider modified 
propagators of more general functional forms. We now act on the exponent on far right with the z derivatives, and 
write z p as id/dp p , integrate by parts over the p variables, and perform the z integral to obtain a delta function over 
p, resulting in: 



9(r) 



t p -m 2 t 

dr 



(27T) 4 



d^ke-ik-ix-v) e -(t-r)u(k 2 ) 



1 d d 



(137) 



i?VA^*ve-™( fe2 V(fc 2 ) 



3 dkP dk x \ pX 

In the second term in the above, we note by the symmetries of the Riemann tensor that the k x derivative can only 

x pv\) 



act on k p (else there will be a term symmetric in k v k contracted with Rg^x), resulting in: 



ft p -rn 2 t 

? (r) = / dr 6 



(2tt)< 



d i kc -tk-(x-y) c -(t-r)u(k 2 ) 9 { U '( k2 ) Rf * k e ~ru(k 2 ) 

dk x \ 3 A M 



which evaluates as: 



g(r) 



(2nY 



d 4 k e~ tu ^ 



u'(k 2 ^ 



>.Jlt + ^tu"(k 2 )R p kpk x - j w ,2 (fc 2 )^fc M fc A 



(138) 



(139) 



where we have performed the trivial r integrals and set x = y. This is a general expression valid for any form of the 
kinetic operator u{k 2 ). We now realize that we can work in a basis in k space where the Ricci tensor is diagonal. 
Furthermore, we note that one can decompose the Ricci tensor in this basis as: 



m = + c p 
11 v 4. 1 



(140) 
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where c£J is diagonal and traceless, and does not renormalize the Einstein Hilbert term (in addition to vanishing upon 
integration in isotropic though not necessarily homogeneous spacetimes). Substituting B% — S^R/4 into the above, 
we then arrive at the expression for a DGP-like propagator: 



Rt 



M 
6 -4fc 



Rt 



1 

12 



Af\ 2 . 



(*+t) 



Rr 



(141) 



from which, it is a straightforward process to compute the term in the heat kernel proportional to the spacetime 
integral of the Ricci scalar as: 



9(r) 



>R 



16tt 2 



l 

&t 



96 



M 
327 



[1 - Erf(MVi/2)Vnie 



M 2 t/4] 



MH 

~6~ 



(142) 



Furthermore, we note that the renormalizations induced by the propagator □ + M\J — □ (i.e. with the other branch 
of the square root) can be similarly calculated, and found to result in the same expression above but with the naive 
substitution M —> —M: 



9{r) 



>R 



16tt 2 



l 

&t 



96 



/M 2 t \ 
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[43] There is also the attitude that discounts the bare cosmological constant problem completely, but attempts to explain dark 
energy and its matter tracking (thus resolving the co-incidence problem) via quitessence, and related models [f| 

[44] Notably, we directly encounter here the fact that diffeomorphism invariance (and the resulting gravitational Ward identities) 
do preclude shifting the zero momentum pole of the graviton propagator, or the generation of extra poles under RG flow 
in the same way that gauge invariance protects gauge boson masses under RG flow via the gauge Ward identities. That 
is, we can consistently generate an effective graviton mass, or extra massive resonances (provided the unitarity bounds on 
a in (O are satisfied) around non-trivial backgrounds under RG flow fully consistent with diffeomorphism invariance 

El- 
See also the footnote at the end of section lib) 

[45] We note that non-local contributions to £ e / / from loop corrections can have several sources, including quantum corrections 
to the quantum state defining the expectation value (|16|) . in addition to contributions which remain finite in the limit x' — > x 
in (|17[) (i.e. are sensitive to the large scale structure of the manifold) [23| . 

[46] Furthermore, we take note here that it is entirely consistent to generate massive spin two poles under RG flow consistent 
with diffeomorphism invariance. This is rather unlike the situation in gauge theories where gauge invariance (via the 
Ward identities) protects gauge boson masses under renormalization. The reason for this is that Einstein-Hilbert gravity 
is a non-renormalizable theory, and corresponds only to the two first terms in a low energy effective theory expansion 
which samples a theory space containing many (in principle infinite) higher dimensional operators under RG flow, unlike 
in renormalizable theories where the theory space consists of a small number of dimensions. Therefore, just as one can 
write down higher order curvature invariants where the graviton propagator can acquire a spin two pole consistent with 
diffeomorphism invariance, so a given RG trajectory flow through regions of theory space where extra poles for the graviton 
propagator result. 

[47] In [26fl. renormalization over a manifold with a boundary results in the leading correction of r 3//2 rather than the r 2 as we 
have in our case. 

[48] The IR cut-off in the proper time integral of the heat kernel has the net effect of suppressing the integration of all 
light momenta. There are several schemes one could have utilized towards this effect, but the hard cut-off remains truest 
to the Wilsonian prescription, and for our purposes has the added advantage of being the easiest scheme to deal with 
calculationally. 

[49] The reader might be under the impression that we measure Gjv via solar system or larger scale observations, but these 
necessarily involve masses of distant objects whose mass we infer from other observations. As the usual gauge theory 
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context, the only way to to experimentally determine the strength of a given interaction at a given scale, and hence fix 
the renormalization conditions, is to scatter particles with known charges off each other and infer the interaction strength. 
The gravitational analog of this would be to test the gravitational force between particles with known masses- i.e. conduct 
Cavendish experiments. 

[50] We also note that unlike the models underlying (|65p . when viewed as a toy example of a free theory, (|65p does not have 

any issues with respect to breaking diffeomorphism or Lorentz invariance. 
[51] Again suspending our disbelief, we consider for a moment the case where the momentum space propagator (|65p can be 

written in the form: (^k — M 2 /kj , with M = 1/0, which describes the propagation of massive plane wave modes in 
Euclidean spacetime- though as before, merely corresponds to a tachyonic instability (m 2 = —M 2 ). In this case, similarly to 
the DGP example one can compute the running of Gjv and A to find: G(fi) = — cfjr{\ ~T- aQ d = A(/v) + V^d^„ 

which would satisfy the requisite properties we require of a degravitating model were it a viable example. 
[52] We note that this theory describes a consistent scalar field which only propagates massless quanta. There is an essential 
singularity in the propagator at zero momentum, which would only be of concern if it implied a negative residue (and would 
therefore indicate the presence of ghosts in our system) . Although the residue at an essential singularity is indeterminate, 

m 2 

one can readily avoid the issue by considering the action in the limit e 2 — > of 5* = J d 4 Xy/—g cj>(e n + ti! □]</>. As we shall 
see, such an action can readily be thought of as encoding the non-perturbative resummations of finite renormalization 
terms. 

[53] In addition, if this non-commutativity is restricted to some hidden (i.e. dark) sector which contains the field <j>, we bypass 

having to worry about Lorentz invariance violating effects at low energies. 
[54] Furthermore, we note with interest that when viewed in momentum space, the non-local action (| T3I) can even be viewed 

as the action for a canonical scalar field modified with a suitable IR cut-off function required to compute its functional 

renormalization group equations [37]]. 
[55] In 40] it was shown that it is not the degree of the differential equation defined by a non-local field theory that determines 

whether or it has a well defined boundary value problem. Rather it is the number of poles, their order, and the residue at 

the poles of defined by the kernel of the kinetic term in momentum space. 
[56] This feature is essentially due to the fact that interactions mediated by massless fields have divergent scattering cross 

sections. 

[57] One might have guessed from dimensional considerations alone that the argument should instead have scaled like M 2 H 2 T 4 , 
where H is some dimensionful parameter characterizing the average expansion of the background. This would follow from 
the fact that in a background with a curvature set by the scale H, we know that G° ~ H 2 f(z). Although this would 
have suited the purposes to which we apply our conclusions to even better, explicit calculations for test sources show that 
factors of H cancel such that indeed, the decay of a homogeneous source to first approximation is independent of the rate 
of expansion of the background (see section 5 of [HJ), with decay of the form (|105p resulting. Perhaps easier seen, is that 
the scaling ~ e~ M H T does not yield a sensible H — ¥ limit in the sense that the degravitation filter clearly still result 
in the decay of a test source (|103f) in a flat background, whereas the ansatz indicates that this would not be the case. 



